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ABSTRACT 

Modern  computational  methods  and  experimental  test  procedures  are  being 
developed  to  determine  more  accurately  the  air-blast  loading  and  dynamic  response  of 
naval  ship  structures  subjected  to  weapon  explosions.  Previously,  a  series  of  air-blast 
trials  on  full-scale  stiffened  panels  typical  of  those  used  in  warship  superstructures  had 
been  conducted  for  conventional  scale  weapons  (under  454  kg  HE )  to  evaluate  new  design 
and  damage -assessment  techniques.  The  present  report  describes  the  experimental  results 
and  the  associated  numerical  study  of  the  panel  response  for  a  longer  duration  loading 
generated  by  the  simulation  of  a  8  kiloton  nuclear  weapon  surface  explosion  called  MINOR 
SCALE.  The  panel,  exposed  to  a  side-on  overpressure  blast  loading  of  345  kPa  peak 
and  200  ms  positive  duration,  experienced  only  slight  permanent  deformation  of  about 
25  mm  with  no  apparent  weld  or  buckling  failure.  The  measured  pressure,  strain, 
acceleration  and  displacement  data  correlated  well  with  the  computed  numerical  results, 

thus  verifying  the  computational  methods.  >  ,  .  >  t 
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1.  INTRODUCTION 

In  the  design  of  naval  ships,  a  capability  to  assess  the  dynamic  response  and  damage 
of  structures  subjected  to  impulsive  loads  from  external  or  internal  air-blasts  and 
underwater  explosions  is  required  not  only  for  the  purpose  of  blast -hardening  the  design, 
but  also  for  accurate  prediction  of  the  mechanism  of  failure  and  the  extent  of  damage 
once  design  thresholds  are  exceeded.  For  these  types  of  explosions,  the  vulnerability  and 
combat  readiness  of  the  complete  ship  structure  must  be  assessed.  The  Canadian  DND 
currently  has  procedures  for  design  of  naval  ship  structures  against  blast  and  shock  loading 
some  of  which  are  based  on  theoretical  and  empirical  approaches  used  for  decades.  The 
Defence  Research  Establishment  Suffield  (DRES)  is  developing  modern  computational 
methods  and  experimental  test  procedures  for  more  accurate  determination  of  blast  and 
shock  loading  and  structural  response  in  order  to  upgrade  and  facilitate  design  analyses. 
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Metal  plates  and  panels  are  fundamental  elements  in  many  military  structures; 
particularly  ships  which  are  composed  primarily  of  stiffened  metal  panels  that  make  up 
the  hull,  decks,  bulkheads  and  superstructure.  The  panels  must  be  designed  to  withstand 
air-blasts  from  either  conventional  or  nuclear  weapon  explosions.  In  regard  to  this 
subject,  DRES  has  been  tasked  by  the  Director  Maritime  Engineering  and  Maintenance 
(DMEM)  and  Director  Maritime  Engineering  Support  (DMES)  to  conduct  a  series  of 
air-blast  trials  on  full-scale  stiffened  panels  typical  of  those  used  in  warship  superstructures 
in  order  to  evaluate  new  design  and  damage-assessment  techniques  being  developed  as 
part  of  the  Shock  and  Blast  research  program  at  DRES  (Project  Number  01  ISA  — 
“Methods  for  Warship  Design  to  Withstand  Blast  and  Shock”).  This  research 
complements  and  extends  DRES’s  previous  experimental  and  analytical  studies  on  the 
effects  of  air  blast  on  masts  and  antennas,  which  were  conducted  in  the  1960’s  to  mid  1970’s. 

The  panel  tests  at  the  DRES  Height -of-Burst  ( HOB)  experimental  range  facility 
were  limited  to  blasts  from  conventional  scale  weapons  (under  454  kg  HE)  giving  loads 
in  the  impulsive  realm  with  relatively  short  time  duration.  The  experimental  data  from 
a  series  of  tests  on  square  plates,  and  the  correlation  with  numerical  results  from  the  finite 
element  code  ADINA  have  been  presented  in  References  1  -  6,  and  a  comprehensive 
summary  of  the  experimental  results  from  tests  on  a  stiffened  panel  is  given  in  Reference  7. 

The  present  report  describes  a  study  of  the  response  of  stiffened  panel  structures 
to  longer  duration  loading  typical  of  that  generated  by  nuclear  weapon  explosions.  On 
27  June  1985,  the  U.S.  Defense  Nuclear  Agency  (DNA)  conducted  a  ground  burst 
simulation  of  an  8  kiloton  nuclear  explosion  called  MINOR  SCALE  at  their  Permanent 
High  Explosion  Test  Site  (PHETS),  White  Sands  Missile  Range  (WSMR),  New 
Mexico.  For  that  test,  an  experiment  was  fielded  to  measure  the  response  and  damage 
of  a  full-scale  stiffened  panel  subjected  to  a  blast  loading  of  345  kPa  (50  psi)  peak 
overpressure  and  200  ms  positive  duration.  The  panel  survived  the  blast  loading, 
experiencing  only  slight  permanent  deformation  of  about  25  mm.  The  test  results  and 
preliminary  numerical  correlations  were  reported  at  the  MINOR  SCALE  symposium  held 
in  February  1986  [  8  ].  The  present  report  contains  a  comprehensive  presentation  of  the 
experimental  data,  a  description  of  the  computational  methods,  and  correlation  of  the 
test  data  and  numerical  results. 
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2.  OBJECTIVES 

To  fill  the  technology  gap  on  the  response  of  ship  stiffened  panel  structures  to  long 
duration  air  blasts,  a  panel  was  fielded  at  the  8-kt  nuclear  ground-burst  simulation  MINOR 
SCALE.  The  project  was  undertaken  in  response  to  a  requirement  by  DMES  under  Task 
DMES  26.  The  objectives  of  the  project  were  as  follows: 

1.  Conduct  measurements  of  the  blast  loading  and  structural  response 
on  a  full-scale  ship's  stiffened  panel  subjected  to  a  long  duration, 
air-blast  overpressure  pulse  of  sufficient  strength  to  produce 
permanent  plastic  deformation. 

2.  Establish  the  pattern  and  mechanism  of  failure  and  the  loading 
function  necessary  to  cause  failure. 

3.  Correlate  the  loading  and  response  measurements  with  design  analysis 
computational  methods. 

Thus,  with  results  and  insight  gained  from  the  experimental  data  and  application  cf  the 
theoretical  and  numerical  analyses,  the  overall  objective  is  to  develop  and  validate  new 
design  and  damage-assessment  techniques  for  ship  panel  structures  subjected  to  long 
duration  air-blast  loading  typical  of  nuclear  weapon  explosions. 

3.  EXPERIMENT  DESCRIPTION 

At  event  MINOR  SCALE,  4.8  kilotons  of  Ammonium  Nitrate-Fuel  Oil  ( ANFO ) 
explosive  material  in  a  hemispherical  fibreglass  container  was  detonated  to  generate  the 
shock  overpressure  loading  equivalent  to  an  8  kiloton  nuclear  explosion.  The  full-scale 
stiffened  steel  panel  fielded  at  MINOR  SCALE  was  identical  in  structural  design  to  the 
panels  tested  at  DRES  (see  Figure  1 ).  The  effective  size  of  the  Tee- beam  stiffened  panel 
was  2.44  m  x  4.57  m  with  6.35  mm  plating  (see  Figure  2 ).  The  panel  was  mounted  flush 
with  the  ground  surface  at  a  radial  location  278  m  from  ground  zero  to  provide  a  blast 
loading  of  345  kPa  peak  side-on  overpressure  and  200  ms  positive  duration.  With 
203  x  203  x  50  mm  sections  welded  to  the  152  x  152  mm  box-beam  perimeter  structure, 
the  panel  was  embedded  in  a  reinforced  concrete  foundation  to  simulate  rigidly  clamp. d 
boundary  conditions.  Photographs  of  the  foundation  and  panel  site  during  and  after 
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the  panel  installation  are  shown  in  Figure  3.  Access  to  the  underside  of  the  panel  via 
a  hatch  cover  and  crawl  space  was  provided  to  allow  for  attachment  of  transducer  sensors 
and  hook-up  of  cables. 

As  shown  in  Figure  2,  the  panel  was  instrumented  at  various  locations  with 
transducers  for  measuring  the  side-on  overpressure,  the  resulting  vertical  displacement 
and  acceleration  of  the  panel,  and  the  strain  on  the  panel  surfaces.  The  strain  was 
monitored  40  mm  from  the  beam  and  outer  boundary  edges.  In  addition  to  the  panel 
instrumentation,  two  DRES  blast  gauges  [9,  10]  were  located  at  either  end  of  the  panel 
(  Figure  3(b))  to  measure  the  free-field  blast-wave  characteristics :  static  and  stagnation 
overpressure,  shock  speed,  and  density.  The  transducer  analog  signals  from  the  panel 
were  cabled  to  an  instrumentation  bunker,  and  recorded  on  magnetic  tape  for  later 
digitization  and  data  processing.  Following  the  blast,  the  panel  was  examined  and 
surveyed  to  determine  the  shape  of  the  deformed  surface  and  the  extent  of  the  damage. 

4.  COMPUTATIONAL  METHODS 

An  important  objective  of  the  project  was  to  correlate  the  blast  loading  and 
structural  response  measurements  with  numerical  predictions.  This  section  will  describe 
the  various  computational  methods  that  have  been  developed  or  are  available  at  DRES. 

4.1  Blast  Loading 

For  the  blast  loading  at  MINOR  SCALE,  the  location  and  orientation  of  the  panel 
flush  with  the  ground  reduced  the  loading  to  a  simple  side-on  overpressure  of  345  kPa 
peak  and  200  ms  positive  duration.  The  loading  was  essentially  spatially  uniform  over 
the  panel  surface  due  to  the  long  duration  of  the  blast  wave  except  for  the  initial  passage 
of  the  shock  front  across  the  panel. 

There  are  various  methods  available  for  describing  the  time  history  of  the  side-on 
overpressure.  Baker  et  al.  [11,  12]  provide  simplified  methods  for  representing  the 
pressure  time-history.  For  the  positive  phase  only,  a  simple  exponential  equation  of  the 
form 


P(t)  - 


(1) 
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where  P,  is  the  peak  overpressure,  b„  the  decay  constant,  and  T,,  the  positive  duration 
is  generally  adequate.  As  an  improvement,  the  modified  Fricdlander  equation 

P(t)  =  P,(l  t  T, )  e 1 '  V  (2) 

has  been  developed  [  13  ].  This  functional  form  gives  a  reasonable  fit  to  experimental 
data  for  the  positive  phase  and  the  negative  phase  up  to  a  time  t  =  1.5  T,,. 

For  structural  response  analysis  using  finite  element  codes,  it  is  advantageous  to 
have  a  computerized  procedure  for  automatically  generating  pressure  loads  resulting  from 
a  specified  explosion.  MARTEC  Ltd.  computerized  Glasstone’s  time  functional 
descriptions  [  14]  of  air-blast  loading  on  simple  2-D  and  3-D  structures,  including  flat 
panels,  vertical  walls,  semicylindrical  arches  and  closed  ‘box- like’  objects.  Tabulated 
values  of  the  dimensionless  air-blast  parameters  (except  for  the  decay  constant  b)  in 
Equation  (2)  were  obtained  from  Reference  15.  The  data  for  the  decay  constant  were 
extracted  from  Reference  11.  The  air-blast  loading  program  developed  by  MARTEC 
is  described  in  Reference  16. 

More  sophisticated  computerized  methods  that  have  been  developed  for  blast  loading 
predictions  include  the  UV1C  Air  Blast  Compendium  [  17  ],  the  Flux  Corrected  Transport 
(FCT)  Method  [18]  and  the  HULL  code  [19],  The  Compendium  generates  time 
functional  loading  from  curve-fits  to  experimental  blast  data  using  the  Random  Choice 
Method  ( RCM )  [20].  The  FCT  and  HULL  codes  determine  the  blast  loading  by  solving 
the  theoretical  equations  describing  the  unsteady  compressible  inviscid  fluid  flow. 

Results  from  these  various  methods  will  be  presented  in  the  next  section. 
Comparisons  with  the  experimental  results  are  given  in  Section  7. 

4.2  Structural  Response 

Numerous  analytical  and  numerical  methods  have  been  developed  to  analyze  the 
static  and  dynamic  response  of  structures  subject  to  various  types  of  loadings.  Various 
finite  difference  and  finite  element  algorithms  have  been  coded  for  predicting  the  static 
and  dynamic,  elasto- plastic,  large  deflection  behaviour  of  beams,  plates  and  stiffened 
panels  [21,  22],  For  linear  and  nonlinear  response  analyses,  the  finite  clement  codes 


UNCLASS1I  IFD 


-•  -C  v‘ 


v 


UNCLASSIFIED 


VAST  [23  ]  and  AD1NA  [24]  have  been  used  at  DRES.  The  code  VAST,  a  linear  code 
developed  by  M ARTEC  Ltd.,  has  an  extensive  finite  element  library,  pre-  and  post¬ 
processors  for  modelling  and  plotting  purposes  [25],  and  a  multi-level 
substructuring /superelement  dynamic  analysis  capability.  Special  programs  have  been 
developed  to  automatically  generate  the  finite  element  model  of  a  stiffened  panel  structure 
[  26  ]  and  to  describe  the  time  history  of  the  pressure  loading  [16].  Translators  [27,  28  ] 
have  also  been  written  to  interface  ADINA  with  the  suite  of  VAST  programs  for  an 
enhanced  analysis  capability  [29].  Gross  deformations  involving  both  material  and 
geometric  nonlinearities  can  be  handled  by  ADINA. 

The  basic  approach  in  any  nonlinear  structural  analysis  is  to  initially  conduct  linear 
analyses  to  establish  a  suitable  finite  element  model  representation  of  the  structure  and 
determine  the  elastic  response  behaviour.  Various  finite  element  models  can  be  generated 
automatically  using  the  VAST  pre-processing  program.  Static  analyses  can  then  be 
performed  to  refine  the  mesh.  Natural  frequencies  and  corresponding  mode  shapes  can 
be  calculated  using  either  direct  or  subspace  iteration.  Dynamic  elastic  response  analyses 
is  generally  performed  using  direct  integration  rather  than  modal  superposition  because 
of  the  impulsive  nature  of  air-blast  shock  loading.  Based  on  the  results  of  the  VAST 
linear  elastic  analysis,  structural  response  can  then  be  analyzed  with  ADINA  to  predict 
nonlinear  static  yielding  and  deformation,  and  the  dynamic  behaviour  leading  to  the  final 
deformation. 

Other  numerical  computational  methods  have  also  been  developed  for  response 
analyses  of  structures  subject  to  blast  loading.  Olson  and  Anderson  at  the  University 
of  British  Columbia  have  developed  a  finite  element  code  FENTAB  [22]  for  beam 
modelling,  and  applied  the  finite  strip  technique  [30]  to  simplify  the  response  analysis 
of  stiffened  and  unstiffened  panel  structures  [31,  32].  These  methods  reduce  the 
computational  requirements  substantially;  thus  making  them  effective  for  preliminary 
design  analyses  and  predictions.  The  finite  element  code  DRESDYN,  developed  by 
Hansen  and  Heppler  at  UTIAS  has  been  tailored  for  response  analysis  of  shell/plate 
structures  exposed  to  a  nuclear  blast  and  thermal  environment  [33,  34],  The  radiant 
thermal  energy  effect  is  represented  in  the  form  of  a  temperature  distribution  initially 
applied  to  the  finite  element  model.  The  response  of  the  structure  subjected  to  the 
pressure-time  history  is  then  analyzed. 
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Structural  response  predictions  using  the  VAST  and  AD1NA  finite  element  programs 
will  be  presented  in  the  next  section.  Final  numerical  results  from  ADINA  as  well  as 
the  UBC  finite  element  and  finite  strip  codes  will  be  compared  with  experimental  data 
in  Section  7. 

5.  PRETEST  DATA 

Prior  to  the  MINOR  SCALE  event,  the  response  and  damage  of  the  stiffened  panel 
were  predicted  based  on  experimental  data  and  analyses  from  the  panel  tests  conducted 
at  DRES  and  finite  element  numerical  calculations.  The  DRES  test  data  were  vital  for 
estimating  the  transducer  response  levels,  understanding  the  possible  modes  and 
mechanisms  of  failure,  and  validating  of  the  finite  element  ADINA  code  for  structural 
analyses  of  the  panel  structures  subjected  to  air-blast  loading. 

5.1  DRES  Panel  Data  and  Analyses 

A  series  of  air-blast  response  trials  of  a  stiffened  panel  (shown  in  Figure  1 )  had 
been  conducted  at  the  DRES  Height -of- Burst  facility.  The  charge  weight  and 
stand-off  height  were  varied,  as  summarized  in  Table  1,  to  investigate  both  the  linear 
elastic  response  and  the  nonlinear  plastic  behaviour  of  the  panel.  The  positive  duration 
of  the  blast-wave  was  generally  less  than  one  third  the  panel  fundamental  period  ( 24  ms) 
so  that  the  loading  was  of  an  impulsive  nature. 

The  first  trial  (No.  311)  was  a  “twang”  test  to  check  procedures  and 
instrumentation.  The  next  two  trials  (No.  312  and  313)  provided  elastic  structural 
response  data.  Trial  No.  313  was  just  below  yield  with  maximum  dynamic  strain  levels 
approximately  twice  the  static  yield  value.  The  fourth  trial  (No.  314)  took  the  panel 
just  above  yield  with  measured  strain  levels  of  4000  pr  and  maximum  permanent 
deformation  of  about  25  mm.  To  achieve  an  even  higher  blast  loading,  two  94  kg  charges 
were  detonated  in  the  fifth  trial.  This  final  shot,  which  provided  a  peak  reflected 
overpressure  of  6375  kPa,  produced  substantial  plastic  deformation  as  shown  in 
Figure  4.  Measured  permanent  deformation  at  the  midspan  of  the  stiffeners  was  about 
90  mm  and  the  center  panel  deformed  an  additional  15  mm.  Maximum  dynamic 
deflections  and  strains  were  approximately  125  mm  and  6200  /us,  respectively.  The  large 
deflections  resulted  in  substantial  membrane  stresses  being  developed  in  the 
plating.  However,  no  visible  cracking  or  weld  failure  was  observed. 
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The  fundamental  natural  frequency  was  observed  to  be  about  41  Hz  (period  24  ms) 
corresponding  to  the  natural  resonance  of  a  single  clamped  panel  segment.  At  event 
MINOR  SCALE  the  duration  of  the  shock  wave  was  expected  to  be  about  200  ms,  so 
that  the  response  of  the  panel  would  be  in  the  dynamic  realm.  Based  on  the  D RES  panel 
results,  particularly  trial  number  314,  it  was  anticipated  that  the  MINOR  SCALE  panel 
with  a  peak  overpressure  of  345  kPa  would  experience  initial  plastic  deformation  of  about 
50  mm  followed  by  elastic  resonance  vibrations  and  a  final  permanent  deflection  of  the 
order  of  25  mm.  The  maximum  strains  and  accelerations  were  predicted  to  be  about 
4000  and  2000  G,  respectively. 

Following  the  panel  trials  at  DRES,  various  linear  and  nonlinear  analyses  were 
conducted  to  validate  the  finite  element  codes  VAST  and  ADINA  for  structural  response 
to  the  blast  wave  loading.  These  analyses  required  tensile  tests  at  quasi-static  and  low 
strain  rates  to  obtain  the  mechanical  properties  of  the  panel  material.  The  following 
values  of  static  material  properties  were  measured :  Young’s  Modulus  =  207  GPa 
(30  x  106  psi);  Poisson’s  Ratio  =  0.3;  Yield  Stress  =  310  MPa  (45  ksi);  Tangent 
Modulus  =  1230  MPa  (  178  x  103  psi);  Ultimate  Tensile  Stress  (at  t  =  0.2)  =  565  MPa 
( 82  ksi).  The  mass  density  was  taken  as  7770  kg/m3.  For  an  anticipated  strain  rate 
of  1.2  s'1,  the  dynamic  Yield  Stress  and  Ultimate  Tensile  Stress  values  were  determined 
to  be  375  MPa  (54.4  ksi)  and  620  MPa  (90  ksi),  respectively.  The  VAST  and  ADINA 
structural  results  have  been  correlated  with  the  test  data  and  reported  in  References  1  -  5. 
For  the  short  duration  type  of  blast  loading  at  DRES  and  for  moderate  deformation  of 
the  panel,  good  agreement  was  obtained  between  the  predicted  and  observed  panel 
motion. 

5.2  Preliminary  Predictions 

Based  on  the  above  calculations,  linear  and  nonlinear  analyses  were  conducted  to 
obtain  more  accurate  predictions  of  the  dynamic  response  expected  at  MINOR  SCALE. 
The  results  of  these  calculations  were  used  for  finalizing  the  design  of  the  experiment  and 
for  setting  the  transducer  signal  conditioning  and  amplifier  gain  to  give  optimum  range 
for  the  transducer  sensitivity  and  recording. 

Using  the  finite  element  model  shown  in  Figure  5  for  a  portion  of  the  stiffened 
panel,  nonlinear  static  ADINA  analyses  [  35  ]  were  performed  (see  Figure  6 ).  Compared 
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to  the  linear  results,  the  nonlinear  geometry  and  material  effects  (i.e.,  large  displacement 
and  plasticity)  substantially  alter  the  predicted  displacement.  The  nonlinear  stiffening 
effect  due  to  membrane  stresses  is  apparent  even  at  low  levels  of  loading;  and  at  210  kPa, 
yielding  of  the  material  becomes  significant.  Based  on  the  nonlinear  material  ADINA 
static  analysis,  the  final  deflection  of  the  panel  for  an  overpressure  of  345  kPa  was  then 
estimated  to  be  about  50  mm  with  maximum  displacements  during  the  first  few  cycles 
of  motion  being  as  large  as  100  mm. 

Calculations  were  also  undertaken  to  predict  the  time  history  of  the  dynamic 
response  behaviour  of  the  panel.  The  expected  pressure-time  history  loading  on  the  panel 
predicted  by  the  simplified  method  (  Equation  2  )  of  Baker  et  al.[  11,  12  ]  and  the  M  ARTEC 
computerized  method  [  16]  were  used  for  selecting  the  site  of  the  panel.  More  detailed 
prediction  of  the  expected  pressure  loading  at  MINOR  SCALE  was  provided  by  S- Cubed 
Inc.  using  the  EIULL  code  [  19],  Figure  7  shows  the  predicted  pressure-time  history  at 
a  location  near  the  panel  site.  The  amplitude  was  scaled  to  obtain  a  pressure-time  curve 
with  a  peak  overpressure  of  345  kPa.  The  positive  impulse  was  predicted  to  be 
15.9  MPa-ms. 

Using  the  loading  data  provided  by  S-cubed  as  input  to  ADINA,  the  dynamic- 
response  of  the  panel  was  computed  for  the  first  40  ms.  The  displacement  -time  histories 
predicted  at  the  center  of  the  panel  (node  289,  transducer  D  12)  and  at  the  midspan  of 
the  beam  (node  303,  transducer  D  13  )  are  shown  in  Figure  8.  The  initial  elasto-plastic 
deformation  due  to  arrival  of  the  shock  front  was  predicted  to  be  -  75  mm  for  the  panel 
and  -  50  mm  for  the  beam.  This  was  followed  by  elastic  vibration  at  1 10  Hz  (9  ms 
period)  which  tended  to  approach  final  deflection  values  of  about  -45  mm  for  the  panel 
and  -35  mm  for  the  beam. 

The  acceleration -time  histories  predicted  for  accelerometers  A  2  and  A  3  located 
at  positions  D  12  and  D  13,  respectively,  are  shown  in  Figure  9.  The  initial  negative 
acceleration  was  expected  to  be  about  900  G,  with  a  subsequent  positive  acceleration  about 
1600  G  for  the  panel  and  1200  G  for  the  beam.  The  prediction  for  acceleration  A  7  at 
node  145  was  similar  to  A  2. 

At  the  time  of  the  calculations,  the  time  history  for  the  components  of  stress  rather 
than  strain  was  available  from  the  ADINA  output.  Because  plastic  deformation  of  the 
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panel  was  expected,  the  stress  could  not  be  directly  converted  to  strain.  For  the  purpose 
of  setting-up  the  strain  amplifiers,  the  AD1NA  stress  predictions  and  the  strain 
measurements  from  the  previous  DRES  trials  were  considered  quite  adequate.  The 
dynamic  stress  predictions  were  converted  to  strain  using  Young’s  Modulus  for  the  elastic 
portion  and  Tangent  Modulus  to  account  for  plastic  yielding.  Strain  levels  of  the  order 
of  4000  iue  were  predicted  on  the  upper  surface  due  to  combined  bending  and  membrane 
effects. 

6.  EXPERIMENTAL  RESULTS 

Immediately  after  the  MINOR  SCALE  Event,  preliminary  plots  of  all  transducer 
signals  were  generated  to  determine  what  test  data  were  successful  recorded  and  to  provide 
data  for  the  ‘quick-look’  report.  The  panel  was  examined  to  assess  the  extent  of  damage. 
The  final  deformed  shape  was  measured  using  a  precision  survey  level.  Subsequently, 
the  recorded  time  history  data  for  pressure,  displacement,  acceleration  and  strain  were 
digitized.  The  computer  program  SPADE  [36]  was  used  to  correct  the  signal  baseline, 
apply  the  calibration  factors,  and  remove  any  extraneous  noise  signals.  This  program 
was  also  used  to  integrate  the  acceleration  to  obtain  velocity  and  displacement,  average 
the  pressure-time  history  traces,  and  to  produce  time  history  plots. 

6.1  Blast  Loading 

Figure  10  shows  the  pressure-time  histories  recorded  by  transducers  P  1,  P  2  and 
P  3  located  at  the  panel  surface  and  P  8  at  a  blast  gauge  station  [9]  near  the  panel.  The 
inserts  show  the  initial  50  ms  of  the  pressure  pulses.  These  pulses  have  some  unexpected 
perturbations  which  could  be  attributed  to  ground  surface  roughness  and  boundary  layer 
effects,  local  anomalies  in  the  shock  front,  ground/panel  surface  discontinuity  or  strain 
gauge  cable  obstruction  near  the  pressure  transducers.  The  peak  overpressure  was 
approximately  345  kPa  (  50  psi),  with  a  positive  impulse  of  15.7  MPa-ms  and  a  positive 
duration  of  200  ms.  At  750  ms,  a  secondary  shock  characteristic  of  blast  waves  for  high 
explosives  is  evident. 

As  evident  from  the  inserts  in  Figure  10,  there  are  differences  in  arrival  time  of 
the  shock  front  at  the  various  transducer  locations.  In  this  figure,  zero  time  is  chosen 
as  the  average  of  the  arrival  times  at  the  four  transducers.  The  differences  in  arrival  times 
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is  partially  due  to  the  velocity  of  the  shock  wave  across  the  panel.  The  shock  wave  transit 
time  from  the  front  to  back  edge  of  the  panel  is  about  3  ms.  However,  correlation  of 
the  shock  arrival  times  from  all  transducers  (reported  in  Reference  9)  shows  that  there 
is  evidence  that  the  incident  shock  front  was  distorted  in  shape  at  this  site  location.  Over 
a  span  of  the  panel  site,  the  shock  front  had  an  inflected  curvature  and  was  obliquely 
incident  by  as  much  as  28  degrees  from  the  expected  radial  propagation.  Since  there 
is  no  above-ground  structures  within  100  m  of  this  site,  a  possible  explanation  for  this 
disturbance  is  a  reflection  from  a  road  and  ditch  alongside  the  panel  site.  The  road  had 
a  hard-packed  gravel  surface  causing  slightly  higher  shock  velocity  and  the  ditch  was  about 
0.7  m  deep  by  2  m  wide.  Both  ran  parallel  to  the  blast  radial  line  and  came  within  about 
3  m  from  the  panel  site. 

Although  some  aspects  of  the  blast  flow  arrival  times  appear  irregular,  these 
irregularities  should  not  strongly  affect  the  panel  response.  Thus,  for  structural  analyses 
purposes,  the  arrival  times  of  the  four  signals  were  synchronized  and  the  amplitudes 
averaged  to  generate  the  representative  pressure- time  history  shown  in  Figure  11.  This 
averaged  pressure  curve  was  used  as  the  applied  load  in  the  final  numerical  analyses 
performed  by  MARTEC  (Section  7.2). 

6.2  Structural  Response 

The  displacement -time  histories  measured  by  transducers  D  12  and  D  13  (see 
Figure  2)  for  the  panel  center  and  beam  midspan,  respectively,  are  shown  in  Figure  12. 
The  initial  deflection  is  approximately  -  53  mm  for  the  panel  and  -  30  mm  for  the  beam, 
followed  by  elastic  oscillations  about  a  mean  that  shows  some  recovery  until  the  arrival 
of  the  second  shock.  It  is  important  to  note  that  the  recovery  of  the  mean  deflection 
follows  a  shape  similar  to  the  exponential  like  decay  of  the  blast  load  amplitude  ( Figure  1 1 ). 
This  in  inherently  related  to  the  long  duration  of  the  MINOR  SCALE  blast  wave.  The 
final  (permanent)  deflection  is  approximately  -25  mm  for  the  panel  and  -  16  mm  for 
the  beam;  corresponding  closely  with  the  post-shot  survey  measurements  of  -  25  mm  and 
-  15  mm  for  the  panel  and  beam,  respectively.  A  photograph  of  the  panel  structure  and 
a  diagram  of  the  final  shape  of  the  panel  surface  ( magnified  by  a  factor  of  10 )  are  presented 
in  Figure  13.  Values  for  the  permanent  deflections  at  the  panel  centers  and  beam 
midpoints  are  included  in  this  figure. 


UNCLASSIFIED 


£ 

.  v  ■-  ,%  /■  .VV’V’VWVvV^V  V  v  v  ">■  V  V.V.V  V.V.V '/.V  V  V,v 


UNCLASSIFIED 


/ 12 


The  accelerations  measured  by  the  various  accelerometers  are  presented  in  Figures  14 
and  15.  The  accelerations  A  1  and  A  3  at  the  midspan  of  the  beams  are  very  similar. 
When  the  traces  of  A  2  and  A  3  are  plotted  together  the  increased  motion  of  the  panel 
center  compared  to  the  beam  midspan  can  be  clearly  seen.  Accelerations  recorded  at 
the  other  points  on  the  pane!  are  similar  to  that  at  A  2.  The  combined  plots  in  Figure 
15  indicate  time  delays  in  the  various  traces  due  to  the  transit  time  of  the  shock  front 
across  the  panel  (about  3  ms)  as  well  as  the  additional  time  difference  caused  by  the  local 
oblique  flow  anomaly  in  the  shock  front  previously  discussed. 

By  integrating  the  acceleration  twice,  the  displacement  motion  is  obtained. 
Figure  16  shows  the  resulting  displacement-time  histories  at  various  points  on  the  panel. 
In  the  upper  two  plots  for  the  panel  center  and  beam  midspan,  respectively,  the  integrated 
displacements  are  compared  with  the  traces  measured  directly  with  the  displacement 
transducers.  Excellent  agreement  is  obtained  in  both  cases,  confirming  the  accuracy  of 
the  displacement  obtained  from  double  integration  of  the  acceleration.  Displacements 
at  the  center  of  the  adjacent  panel  segments  given  by  A  4  and  A  5  are  similar  in  level 
to  A  2.  but  A  5  exhibits  slightly  more  elastic  vibration,  being  near  the  panel  edge. 

Plots  of  various  measured  strain-time  histories  are  presented  in  Figures  17  -  19. 
No  data  were  available  from  S  15,  S  17,  S  21  and  S  25  because  of  gauge  failure.  The 
strain  measurements  on  the  upper  and  lower  surfaces  shown  in  Figure  17  illustrate  the 
tensile  and  compressive  bending  strains,  and  the  presence  of  an  in-plane  membrane  strain 
corresponding  to  the  difference  between  the  two  signals.  The  strain  signals  presented 
in  Figures  18  and  19  show  the  consistency  obtained  across  the  panel,  and  the  relative  levels 
of  the  strain  normal  and  tangential  to  the  panel  edge  about  40  mm  from  the  boundary.  In 
the  bottom  diagram  of  Figure  19,  the  initial  signals  from  strain  gauges  S  18  and  S  20  arrive 
at  slightly  different  times  (about  3  ms)  due  to  the  transit  time  of  the  shock  front  across 
the  panel.  Peak  strain  levels  of  the  order  of  3500  pc  were  measured  by  the  gauges  on 
the  upper  surface  due  to  combined  bending  and  membrane  effects.  Along  the  actual 
lines  of  attachment  of  the  plating  to  the  beam  stiffeners  and  outer  boundary  edges,  even 
higher  levels  of  strain  arc  anticipated  due  to  plastic  hinge  effects.  The  initial  strain  rate 
calculated  from  the  plots  ranges  from  0.76  s  '  to  1 .9  s'1  which  spans  the  anticipated  value 
of  1.2  s  '  used  'o  determine  the  material  properties  in  the  pretest  finite  element  analysis 
(see  Section  5.2). 
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7.  NUMERICAL  RESULTS  AND  CORRELATIONS 

7.1  Blast  Loading 

Since  the  dynamic  response  of  the  panel  structure  depends  on  the  loading  function, 
the  measured  pressure-time  history  was  employed  for  the  final  numerical  calculations. 
Although  the  finite  element  code  can  incorporate  the  transit  time  of  the  shock  front  over 
the  panel  and  any  other  spatial  variations,  the  analysis  was  simplified  by  assuming  that 
the  pressure  load  was  applied  instantaneously  and  that  it  was  spatially  uniform  over  the 
panel  surface.  The  pressure-time  history  {shown  in  Figure  11)  was  generated  by 
synchronizing  and  averaging  the  amplitudes  of  the  actual  measured  pressure  traces.  The 
averaged  loading  data  agree  well  with  the  HULL  code  predictions  (  Figure  7  )  except  for 
the  perturbations  in  the  initial  waveform  between  2  to  14  ms.  As  will  be  evident  by  the 
following  results,  the  perturbations  in  the  initial  waveform  do  affect  the  dynamic  response 
and  vibrational  characteristics  of  the  panel  structure. 

7.2  Finite  Element  (ADINA)  Results  (MARTEC) 

Using  the  ADINA  code  and  a  refined  finite  element  model  (Figure  20), 
MARTEC  [  37  ]  re-calculated  the  displacement  and  acceleration  at  the  center  of  the  panel 
(node  457,  transducers  D  12  and  A  2)  and  at  the  midspan  of  the  beam  (node  472, 
transducers  D  13  and  A  3 )  (see  Figures  21  and  22,  respectively).  As  shown  in  Figure  21, 
the  initial  elasto- plastic  deformations  are  now  determined  to  be  -48  mm  for  the  panel 
and  -28  mm  for  the  beam;  about  40%  smaller  than  the  preliminary  predictions  (see 
Figure  8)  based  on  the  HULL  code  predicted  loading.  The  reduction  in  maximum 
dynamic  displacements  is  a  result  of  the  reduced  effective  loading  due  to  the  perturbations 
in  the  actual  shock  wave  between  2-14  ms.  By  projecting  the  curves  to  later  times, 
it  appears  that  the  permanent  deflections  approach  about  -25  mm  for  the  panel  and 
-  15  mm  for  the  beam,  corresponding  closely  with  the  measured  deformations  of  the  panel 
(  -25  mm)  and  beam  (  -  16  mm)  shown  in  Figures  12  and  13.  The  acceleration  time 
history  for  accelerometer  A  2  (Figure  22)  has  a  range  (  -900  G  to  1600  G)  similar  to 
the  preliminary  predictions  (see  Figure  9),  but  the  acceleration  at  A  3  on  the  beam  has 
a  range  (  -500  G  to  800  G)  which  is  somewhat  less  than  that  originally  predicted. 

Further  comparison  of  the  ADINA  numerical  results  with  the  experimental  motion 
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data  shows  good  agreement  (  Figures  23  and  24 ).  The  measured  and  predicted  acceleration 
and  displacement  at  the  panel  center  and  beam  midspan  have  similar  levels  and 
waveform.  The  initial  plastic  deformations  occur  within  the  first  6  ms  of  the  applied 
shock  wave  loading  and  peak  deflections  occur  within  one  period  (24  ms).  The 
experimental  data  also  exhibit  some  higher  frequency  vibrational  characteristics  due  to 
harmonic  vibrations  and/or  transverse  wave  vibrations  in  the  plating.  Some  of  these 
\  ibrations  could  also  be  extraneous  noise  generated  by  transducer  ringing,  and/or  vibration. 

7.3  Finite  Element  (ADINA)  Results  (DRES) 

Correlation  of  the  numerical  predictions  with  the  experimental  data  depends  of 
the  accuracy  of  the  finite  element  model  and  the  loading  function.  The  size  of  the  time 
step  for  numerical  integration  and  the  interpolation  procedures  and  convergence  criteria 
are  also  important. 

Using  the  latest  version  of  ADINA  (  1984),  further  analyses  were  conducted  by 
refining  the  model,  extending  the  numerical  calculations  to  1000  ms  and  producing  strain¬ 
time  history  results  for  direct  correlation  with  the  experimentally  measured  data.  The 
new  finite  element  model  (Figure  25)  was  developed  using  the  4-node  quadrilateral 
plate  shell  element.  An  extended  pressure  curve  obtained  by  averaging  transducers  P  1 
and  P  3  was  used  as  the  pressure  loading  (  Figure  26). 

The  excellent  correlation  of  the  displacement  results  for  the  panel  center  well  beyond 
the  second  shock  is  shown  in  Figure  27.  The  prediction  of  too  large  a  final  beam  deflection 
is  due  to  over  stiffness  of  the  finite  element  model  for  the  beam  after  plastic  deformation 
and  the  lateral  membrane  restraint  for  the  remainder  of  the  panel.  Modelling  one  quarter 
of  the  complete  panel  and  refining  the  model  with  more  elements  through  the  web  height 
and  across  the  flange  would  reduce  the  stil  lness  and  improve  the  displacement  predictions 
during  the  recovery  phase  of  the  panel  deflection. 

The  numerical  results  for  strain-time  history  arc  determined  at  the  gauss  integration 
points  of  the  nearest  finite  elements  and  then  interpolated  to  strain  gauge  locations. 
Correlation  of  the  strain -time  histories  for  a  number  of  the  gauges  are  shown  in 
I  igurcs  28  30.  Numerical  results  for  S  15  and  S  17  are  included  in  Figure  28  even 

tough  no  experimental  data  were  available  because  of  gauge  failure  at  these  locations. 
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The  difference  in  start  time  between  the  predicted  and  measured  strain-time  histories  in 
Figures  29  and  30,  is  simply  due  to  the  sweep  time  of  the  shock  front  ovct  the  panel  relative 
to  the  synchronized  and  averaged  time  of  the  loading  function  used  for  the  analyses. 

7.4  Structural  Response  (UBC) 

Numerical  results  for  the  panel  motions  were  also  obtained  using  newly  developed 
finite  element  and  finite  strip  codes.  The  finite  element  equivalent  beam  and  finite  strip 
models  of  the  stiffened  panel  are  shown  in  Figure  3 1 .  These  calculations  were  performed 
using  the  pressure-time  loading  shown  in  F  igure  26. 

The  numerical  displacement  results  obtained  using  these  models  for  the  initial  50  ms 
are  illustrated  in  Figure  32  together  with  the  experimental  data  and  the  DRES  AD1NA 
calculations.  The  finite  strip  calculations  were  then  carried  out  to  440  ms  (see  Figure  33  ). 
The  panel  center  is  predicted  to  exhibit  greater  elastic  vibrations  and  then  a  buckling 
instability  at  280  ms,  referred  to  as  an  ‘oil-canning’  phenomena,  after  which  the  plating 
vibrates  in  a  bowed-up  shape  relative  to  the  deformed  beams.  If  the  buckling  instability 
is  ignored  (since  it  did  not  occur),  the  finite  strip  results  are  in  good  agreement  with  the 
measured  data.  The  effects  of  damping  on  the  elastic  vibrations  and  the  occurrence  of 
the  ‘oil-canning’  phenomena  warrant  further  investigation.  There  are  obvious  advantages 
to  using  the  simplified  methods  in  that  computational  times  are  between  10  to  20  minutes 
compared  to  30  hours  for  the  ADINA  calculations,  so  that  the  predictions  are  obtained 
at  substantially  reduced  cost. 

8.  CONCLUDING  REMARKS 

Event  MINOR  SCALE  provided  a  unique  opportunity  to  measure  the  response 
of  the  stiffened  panel  structure  for  a  long  duration  pressure  loading  typical  of  a  nuclear 
weapon  explosion.  The  345  kPa  peak  overpressure  level  with  about  200  ms  positive 
duration  exposed  the  panel  to  a  loading  sufficient  to  cause  permanent  deformation. 

The  permanent  deformations  of  25  mm  for  the  panel  and  15  mm  for  the  beam 
were  less  than  the  pretest  predictions.  Acceleration  levels  between  1000  -  2000  G  were 
about  what  were  expected.  Strain  levels  up  to  3500  \jl  were  measured  on  the  upper  surface 
due  to  combined  tensile  bending  and  membrane  effects.  Higher  strains  are  likely  (but 
were  not  measured)  closer  to  the  boundary  as  a  result  of  the  plastic  hinging  effects. 
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The  discrepancy  between  the  experimental  data  and  the  pretest  predictions  was 
directly  related  to  the  loading  function  used  in  the  analyses.  By  using  the  actual  measured 
pressure  loading  the  dynamic  response  analysis  showed  excellent  agreement  between  the 
numerical  results  and  the  measured  acceleration  and  displacement  data.  This  illustrates 
that  in  order  to  examine  the  details  of  structural  response,  it  is  necessary  to  use  the  actual 
pressure  loading.  Nevertheless,  for  the  purpose  of  preliminary  design  analysis  and 
planning  an  experiment,  the  HULL  code  prediction  was  found  to  be  quite  adequate  in 
giving  the  general  features  of  the  pressure  loading. 

The  success  of  the  panel  experiment  at  event  MINOR  SCALE  substantiated  the 
trials  procedures,  data  recording  and  processing  techniques  that  were  used.  The 
correlation  of  results  verified  the  suitability  of  the  finite  element  code  ADINA  for  detailed 
analysis  of  structural  response,  whereas  the  simplified  computational  methods  based  on 
the  finite  element  equivalent  beam  and  finite  strip  techniques  have  been  shown  to  be  suitable 
for  preliminary  design  analysis. 

The  experiment  fulfilled  the  objectives  of  the  task  by  subjecting  a  full-scale  ship's 
stiffened  panel  to  a  long  duration  air-blast  loading  sufficient  to  cause  permanent  plastic 
deformation,  and  by  establishing  the  resulting  pattern  and  mechanism  of  failure.  The 
good  correlation  of  the  experimental  data  and  numerical  results  validated  the  computational 
methods  for  long  duration  type  of  loading  function  and  for  structures  sustaining  limited 
amount  of  permanent  deformation.  The  methods  can  now  be  used  with  confidence  in 
•'hip  design  for  carrying-out  structural  response  analyses  and  damage  assessment  studies. 

However,  tor  structures  experiencing  considerably  higher  plastic  deformation  and 
damage,  further  te->ts  and  analyse'  are  recommended  for  complete  validation  of  the 
computational  methods. 
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